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Thomas SotiriotQ and Enrico Barauss^j] 

SISSA, International School for Advanced Studies, Via Beirut 2, 34014 Trieste, Italy and INFN, Sezione di Trieste 

(Dated: February 7, 2008) 

The Parametrized Post-Newtonian expansion of gravitational theories with a scalar field coupled to the Gauss- 
Bonnet invariant is performed and confrontation of such theories with Solar system experiments is discussed. 

PACS numbers: 04.25.Nx, 04.50.+h, 04.80.Cc 



■ 

o ■ 
o : 

<N . 

5h ' 

<: 
o ■ 
<n : 

(N '• 
> ■ 

in ; 

\o . 
o ■ 

(N ■ 

\o : 
o ■ 

"o ! 

cr 

i ' 



X 



The observed late time accelerated expansion of the uni- 
verse combined with the problems that arise when trying 
to rationalize it with the simplest of explanations, such as a 
cosmological constant yfl, have triggered an increased inter- 
est in finding alternatives for the nature of dark energy. Scalar 
fields, widely used in the inflationary paradigm [3], constitute 
a familiar way of providing accelerated expansion. Scalar- 
tensor theory is therefore an appealing candidate as an alter- 
native theory of gravity that can provide the desired cosmo- 
logical dynamics Jj]. However, there are motivations from 
String theory to believe that scalar fields might not be cou- 
pled to the Ricci scalar, as in scalar-tensor theory, but to the 
Gauss-Bonnet term, 

g = R 2 -AR^R^ + R^R^x, (i) 

which is topologically invariant in four dimensions. 

To be more precise, one expects to find two types of scalar 
fields in the low energy effective action of gravity coming 
from heterotic String theory: moduli, 0, which are related 
to the size and shape of the internal compactification mani- 
fold, and the dilaton a, which plays the role of the string loop 
expansion parameter. There are reasons to believe that mod- 
uli generally couple to curvature square terms [5] but mod- 
uli dependent higher loop contributions, such as terms cubic 
or higher order in the Riemann tensor, vanish, leaving a cou- 
pling with a Gauss-Bonnet term to be of specific interest J5[0] . 
On the other hand the dilaton usually couples to the the Ricci 
scalar and consequently to matter in the Einstein frame, and 
there are reasons to believe that it might evolve in such a way 
as to settle to a constant J3]. Under these assumptions the 
effective low energy gravitational action takes the form 

JL-^d^-v^+mg , 
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(2) 

where k 2 — 8irG, A is +1 for a canonical scalar and —1 for a 
phantom field (c = % = 1). A straightforward generalization 
of the action is to include a kinetic term and a coupling with 
g for the dilaton a. This will not concern us here but we will 
discuss how we expect it to affect our results. 

Remarkably, it has been shown that action (01 can lead to a 
theory of gravity with desirable phenomenology, including in- 
flation and accelerated expansion 08||9|]. Such a theory seems 
to fit observational data related to cosmology adequately [10] 
and avoid past and future singularities |6|,|8j]. However, a grav- 
itational theory which can pose a viable alternative to General 
Relativity should also have the correct Newtonian and post- 
Newtonian limits, since Solar system tests provide stringent 



constraints. This is the issue that will concern us here. We 
will not consider the case where A = and the scalar field has 
no kinetic term in the action. Such actions are dynamically 
equivalent to an action with a general function of Q added to 
the Ricci scalar and their Newtonian limit has already been 
considered in 111 ill . 

Let us start by reviewing the field equations one derives 
from action (01 1 8]. For the metric we have 
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^-r^ -T-2(V„V„/(0))i2- 
-2g^{Uf{<f>))R - 4(VV M /(0)Jk, - 
-4(V p V u f(<t>))R pp + 4(D/(0))iV + 
+4g^(V^f(^))R pa - 4(VV CT /(^))i? p 

where G^ = R^ - \Rg pv and 

= A (-dpWvt - \g» v d p 4>d P 4\ - \gp»V{4), (4) 

and for the scalar field 

xa < p-v'( ( f > ) + f'(i > )g = o, 



(5) 



where A'(x) = dA/dx and □ = .g^V^V^. 

Let us bring eq. © to a more suitable form for our pur- 
poses. Taking the trace of eq. (01 one gets 

R = k 2 [-T -T* + 2(Of(<p))R - 4(VV < 7(^))i? p(T ] , 

(6) 

where T = g^T pu and = g^T*,. Replacing eq. © 
back in eq. (01, the latter becomes: 
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+2(V„V„/(0))i2 - g^{Vf{(t>))R - 
-4(V%/(^))^ - 4(VV„/fa))ii /v 
+4(n/(0))iV + 2g^ p V a f{<j } ))R pa - 

-4(VV ct /(0))^p- 



(7) 



Following [12 1 ] we can choose a system of coordinates 
in which the metric can be perturbatively expanded around 
Minkowski spacetime. Therefore we write the metric as 
gfiu = Vfiu + h pl/ and the scalar field as <fi = <fio + 5<fi, where 
the value of </>o is determined by the cosmological solution. 
The perturbed field equations are 
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(50)5$ - V'bMW ~ 2^"'(0o) (5cf>) 2 + f'((f> )G = 0(8ct?,8(t>{h^f, h^fo, h^fa), (8) 



^A 4>l - ^V(<p ) + \v\<M5<t> (-1 + /ioo) + 



A[nflat<5' 
R 2 |roo 

+f'{M [2(d d 54>)R + (Dflat^i? - 8(dPd 5cf>)R 0p + 4(D flat ^) J Roo - 2(c^ ff <50)i? p(T 

-4(G>"a ff ^)i? 0p o CT ] } + <W 2 V. W\j>o54>, 0o V, V(^o)/ioo), 
k 2 T 0i + O^Stphfw, 5(f) 2 ,Th 0i , foSfa foli^, V((f> )h 0l ), 

+ 0(5(f>h^ v ,8(f> ,Thij,4>oh^ u , V(<po)fkj), 
I 



r« + o% (- T + ^'(^o)¥ + ^(0o)) 



(9) 
(10) 

(11) 



where Dfl a t denotes the D' Alembertian of flat spacetime. No- 
tice that, as usually done in scalar-tensor theory lfl3i , we 
have neglected all of the terms containing derivatives of (f>o 
multiplying perturbed quantities (e.g. <fio5<f>). This is due to 
the fact that (f>o changes on cosmological timescales and con- 
sequently one expects that it remains practically constant dur- 
ing local experiments. Therefore its time derivatives can be 
neglected as far as solar system tests are concerned. This can 
easily be verified by some order-of-magnitude analysis. Take 
for instance Eq. (|9): the terms containing a time derivative of 
4>q multiplying a perturbation are 0(f((f>o)h fJllJ /(r 2 M 2 )) and 

O(<po8(j>/M%), where O - H M p and / ~ H§ (M p = 1/k 
is the Planck mass and Hq the present Hubble constant) and 



'oo 



rScj) ~ hoi/v 



2 5(j)/v ~ GMq/t (r is 
= y/GM Q /r). On the 



the distance from the Sun and v 

other hand, the C(w 4 ) post-Newtonian correction to Rqq is 
~ (GMq) 2 /r 4 ~ lO- 55 O(/(0 o )^/(r 2 A/ p 2 ),0 o ^/M p 2 ) 
even if r is taken as large as 1000 AU. Therefore, the cor- 
rections coming from terms containing time derivatives of </>o 
multiplying perturbations are at least 55 orders of magnitude 
smaller than post-Newtonian corrections, and neglecting these 
terms cannot affect our results in any way. A similar treat- 
ment applies to the terms containing the potential V multi- 
plying perturbed quantities (s.g.V(4>o)hoo): in order to give 
a reasonable cosmology, V(<fio) should be of the same order 
as the energy density of the cosmological constant and these 
terms cannot therefore lead to any observable deviations at 
Solar system scales. 

We will instead adopt a different treatment required for the 
simple V(<j>o), \(f> 2 terms appearing in the field equations: as 
they need to be of the same order as the energy density of the 
cosmological constant, they will not lead to any observational 
consequences (see lfl5ll and references therein). However, for 
the sake of the argument we will keep track of them: due to 
their small value we can treat them as 0(v 4 ) quantities follow- 
ing lfl5tl and so they will not appear in the 0(v 2 ) equations. 
As far as terms related to V'(4>o) are concerned, we intend to 
just keep track of them for the time being and discuss their 
contribution later on. 

If we now expand the perturbations in the metric and the 



scalar field in post-Newtonian orders, keeping in mind that the 
Parametrized Post-Newtonian expansion requires <f> and /loo to 
0(v 4 ), to 0{v 2 ) and h oi to 0(v 3 ), we can write 
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25(p + 4,dq> . . 
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where the subscript denotes the order in the velocity, i.e. quan- 
tities with a subscript 2 are 0(v 2 ), quantities with a subscript 
3 are 0(v 3 ), etc. So to order 0(v 2 ) this gives 



AV 2 ( 2 <50)-^"(0 O )2^ = O 



(16) 



where V 2 = Sijdidj. Note that since the metric is flat in the 
background Q = (D(v A ). We want <j> to take its cosmological 
value at distances far away from the sources. This is equiva- 
lent to saying that the perturbations due to the matter present 
in the Solar system should vanish at cosmological distances, 
and this can be achieved by imposing asymptotic flatness for 
the solution of eq. (TToT l. i.e. 25(f) — > for r — ► 00. This implies 
25<i> = 0. 

To order 0(v 2 ) for the components 00 and ij and 0(v 3 ) 
for the components Qi, and after applying the standard gauge 
conditions 
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fan ,, ~ h^, n — ~hn 



(17) 



the field equations for the metric take the form 

-V 2 ( 2 /ioo) 

-v 2 ( 2 M 

^ ^V 2 ( 3 /i 0l ) + ~ (2/100 jo) 

which remarkably are exactly the same as in General Relativ- 
ity ni2ll . The well-known solutions are 
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where following [12] we define 

U = G f d 3 x' 



3 , p{x',t) 



\x — x'\ ' 

r: = a / ^K^My,*) 

|sc — ar | 



(24) 
(25) 



w ^ /" , 3 tyvH^^c - x') k (x - x')j 
W t = Gjdx j— (26) 

We already see that the theory has no deviation from Gen- 
eral Relativity at order 0(v 3 ): in particular it gives the cor- 
rect Newtonian limit. It is now easy to go one step further and 
write down the perturbed equations that we need to 0(v A ). 
For the scalar field, using 2 5<fi = 0, we get 



AV 2 ( 4 c 



V"{4> ) 4 6<f> + iQ = 



(27) 



with 



aQ = (2^-00, y)(2^00,ij) — (2^00,ii)(2^00,ij) + {ihijjj) + 
+ {2hij > kl)(2hij ! kl) — (2^ij : fefe)(2^ii,fcfc) — 

—2( 2 hij : ki){2hiijk) + \ihij t u)\2h,kl,ij) , (28) 

where we have again applied the gauge conditions (T% . Using 
eqs. (fJTJ and d22l . eq. (f28b becomes 

4 £ = 8[/ fci £/ fci -8(£/ fcfc ) 2 . (29) 

The solution of eq. d27| i is therefore 



47T 



\x — x'\ 



(30) 



The time-time component of the perturbed field equations 
for the metric to 0{v A ) is 
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(4T00) + -UT) - ~{ 2 h 00 ){ 2 T) 
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1, 



V(&,)( 4 <ty) - 7^(00) 
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(31) 



where we have already used the fact that 2 S<p = 0. Note also 
that no contribution coming from the coupling between <j> and 
the curvature terms in eq. (0 is present in the above equa- 
tions. This should have been expected since in eq. (0 these 
terms always have the structure of two derivatives of times a 
curvature term, and so, due to the fact that in the background 
the metric is flat and </>o is slowly varying, they can only con- 
tribute to orders higher than 0(v A ). 

Let us discuss the contibution of the term proportional to 
V'((/>o). Using eqs. d30b and d29l we can write this term as 
an integral over the sources times a dimensionless coefficient 
k 2 V (0o)/'(0o)- O ne can argue the V'(<f>) should be prac- 
tically zero as far as the post-Newtonian expansion is con- 
cerned ifLil a. This is equivalent to saying that the cosmo- 
logical solution corresponds to a minimum of the potential. 
Even though such assumptions are not exact, they are accu- 
rate enough for our purposes. Note that even in cases where 



V does not have a minimum, well motivated models usually 
introduce exponential forms for the potential and the coupling 
function, i.e. V = V e~ aK ^ and / = /oe bK< ^ where a, b and 
/o are of order unity whereas Vo is as small as the energy den- 
sity of the cosmological constant in order to guarantee that 
the theory will fit observations related to the late time cos- 
mological expansion. This implies that, since k 2 ~ 1/M 2 , 
k 2 V(0o)/'(0o) is dimensionless and of the order of the now 
renowned 10~ 123 . Therefore, we will not take the term pro- 
portional to V'(ipo) into account for what comes next but we 
will return to this discussion shortly. 

We can use the solutions for 2 hoo and 2 hij, the gauge con- 
ditions ( fTTI i and the standard post-Newtonian parametrization 
for matter 111211 to write eq. (l3~TT > as 



V 2 ( 4 /i o + 2C7 2 -8* 2 ) = 



2p ( v 2 - U + -n 



3p 
2p 



V(<po) + 2 Ml 



(32) 



where IT is the specific energy density (ratio of energy density 
to rest-mass density) [12] and 



$2 = G / d 



l3 ^ p(x',t)U(x',t) 
\x — x'\ 



The solution to this equation is 

4 h 0Q = 2U 2 + 4$x + 4$ 2 + 2$ 3 + 6$ 4 



(33) 



where 
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3 , p(x',t)v(x',t) 2 
( ' ! dx • 
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j3 ^ p{x',t)U(x',t)v(x',t) 2 
I x — x I 



Gld s x' Pix '^ 



(34) 

(35) 
(36) 
(37) 



\x — x'\ 

Therefore the metric, expanded in post-Newtonian orders is 

-1 + 2U- 2U 2 + 4$i + 4$ 2 + 



.9oo 



-2$ 3 + 6$ 4 + !L. ( V (ct>o) ~ ) -if. (38) 



7 1 

90j = ~2 Vi ~2 Wi ' 
= (l + 2E/)<%, 



(39) 
(40) 



which, apart from the term related to V (<f>o) — 1/2A0 2 ,, is ex- 
actly the result that one obtains for General Relativity. This 
term corresponds to the standard correction normally arising 
from a cosmological constant, and since V(4>o) — 1/2A0 2 , 
should indeed be of the same order as the energy density of 
the cosmological constant, the contribution of this term is neg- 
ligible on Solar system scales. Since the metric is written in 



4 



the standard PPN gauge one can read off the PPN parame- 
ters 11211 . The only non-vanishing ones are 7 and j3, which 
are equal to 1. Therefore, the theory discussed here seems 
to be indistinguishable from General Relativity at the post- 
Newtonian order. 

The above implies that a gravitational theory with a scalar 
field coupled to the Gauss-Bonnet invariant trivially satisfies 
the constraints imposed on the post-Newtonian parameters by 
Solar system tests. This appears to be due to the fact that the 
terms arising in the field equation for the metric from the cou- 
pling between the scalar and Q in the action always have the 
structure of two derivatives of / times a curvature term. Such 
terms do not contribute to the post-Newtonian expansion to 
0{v A ). This is not the case for other possible couplings of a 
scalar to a quadratic curvature term, such as <fiR 2 . Remark- 
ably, the characteristic structure of such terms can be traced 
back to the special nature of Q, i.e. to the fact that it is a topo- 
logical invariant in four dimensions. 

However, notice that our result strongly depends on the as- 
sumption that V(<f>o) and V'(4>o) are reasonably small so as 
to give a negligible contribution. This assumption stems from 
the fact that V((j>o) will play the role of an effective cosmo- 
logical constant if the theory is to account for the late time 
accelerated expansion of the universe and should therefore be 
of the relevant order of magnitude. Additionally we expect 
that V'(4>q) will also be small enough so that its contribution 
can be considered negligible, based on the fact that either the 
field approaches a minimum at late times, or that the poten- 
tial is of the form V = Voe^ aK ^, where a is of order unity, 
and therefore V (<f>o) ~ kV(4>q), which seem to be true in all 
reasonable models. 

Another important aspect that has to be stressed is that the 
value of f'(cj>o) or f"(4>o) should be suitable so that the post- 
Newtonian expansion remains trustworthy. From eq. ( f30b we 
see that non-trivial corrections will indeed be present in post- 
post-Newtonian orders and, even though such corrections are 
normally subdominant, if f'((f>o) or f"{4>o) are sufficiently 
large they can become crucial for the viability of the theory. 
This was first observed in iflfjll where the same theory, but 
without a potential V, was confronted with Solar system ob- 
servations, considering a nearly Schwarzchild metric as an ap- 
proximation. As mentioned before, the potential plays the role 
of an effective cosmological constant if one wants a theory 
that leads to a late time accelerated expansion as in Q8|, [9[] . If 
this potential is not present, it is the coupling /(</>) between 
the scalar field and the Gauss-Bonnet term that will have to 
account for this phenomenology. In this case it turns out that 
f"(4>o) has to be of the order of the inverse of the cosmolog- 
ical constant, and this is enough to make the post-Newtonian 
approximation break down. Fortunately, models with a poten- 
tial do not suffer from this problem, and as a matter of fact / 
is usually assumed to be of the form / = f^e bK ^ where both 
/o and b are of order unity. So, as predicted also in 11611 . all 
reasonable models with a potential will pass the Solar system 
tests. 

Finally, let us discuss the possibility of including a second 



scalar field in the action, coupled to the Gauss-Bonnet invari- 
ant, which could, for example, be the dilaton. If this second 
scalar is not coupled to matter or to the Ricci scalar, then it 
can be treated using the same approach. If the coupling func- 
tion with the Gauss-Bonnet invariant as well as the potential, 
if present, have similar properties as those discussed above 
we expect our result to remain unaffected. Of course there is 
also the possibility that the dilaton is coupled to matter. This 
goes beyond the scope of the present paper since in this case 
the theory would be phenomenologically different not only in 
what regards Solar system tests but also in other aspect such 
as cosmology, covariant conservation of matter, equivalence 
principle (e.g. see ref. 111711 ) etc. 
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